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Charge ratio of muons from atmospheric neutrinos 
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We calculate the intensities and angular distributions of positive and negative muons produced by atmospheric 
neutrinos. We comment on some sources of uncertainty in the charge ratio. We also draw attention to a poten- 
tially interesting signature of neutrino oscillations in the muon charge ratio, and we discuss the prospects for its 
observation (which are not quite within the reach of currently planned magnetized detectors). 



The MINOS detector under construction in 
the Soudan Mine for study of neutrino oscillations 
with a neutrino beam from Fermilab will be able 
to measure the momentum and charge of muons 
in the energy range from just below one GeV to 
about 70 GeV/c. In addition to its use as the 
far detector for accelerator neutrinos, MINOS is 
sufficiently large to detect a significant number 
of atmospheric neutrinos. Muons entering from 
outside the detector as well as muons originating 
from vertices inside the detector will occur. These 
events will provide a useful calibration beam for 
the detector. 

The detector, which is now under construc- 
tion at Soudan, has already detected atmospheric 
muons and neutrinos. After completion, MINOS 
will continue to collect cosmic-ray data for some 
time before the accelerator neutrino beam turns 
on and afterwards as well. Because the atmo- 
spheric neutrino beam has comparable fluxes of 
neutrinos and anti-neutrinos, it will be possi- 
ble [^y^ to search for any sign of differences in 
oscillation properties of and Pfj,. Such differ- 
ences could arise in theories in which there is 
intrinsic violation of CPT [||. They could also 
arise from matter effects, for example with mix- 
ing of three active neutrino flavors if sin^ 9i2 is 
large enough ||^ or with non-maximal ^ 
mixing. Q 

The purpose of this letter is to calculate the 
rates of atmospheric neutrino interactions in a 
detector like MINOS as a function of direction 
and energy separately for neutrinos and anti- 



neutrinos. We briefly discuss some sources of sys- 
tematic uncertainty, and we note a potentially in- 
teresting effect of oscillations as reflected in the 
charge ratio of atmospheric muons. 

To illustrate the basic result, we start with 
the atmospheric neutrino spectrum of Ref. [Q. 
We integrate the spectra of atmospheric neu- 
trinos folded with da/dy to obtain the rate of 
events with contained vertices. To obtain the 
rates of external muons, the convolution also in- 
cludes the muon range. For simplicity we 
use the GRV94 structure functions |10 to ob- 
tain the spectra of neutrino- induced fj/ and ji^ . 
(We have checked that the results for the lim- 
ited range of x and relevant here are es- 
sentially the same as those obtained with more 
complicated structure functions.) In Table 1 we 
show expected event rates as a function of muon 
energy with and without oscillations, separately 
for contained vertices and for external upward- 
moving neutrino-induced muons. Full mixing 
with Am^ = 2.5 x lO^'^ eV^ ||ll| was assumed for 
the oscillation case. Contained vertices are shown 
integrated over all directions; neutrino-induced 
muons are summed over the upward 27r steradian. 

The first point to note is that the charge ratio 
for neutrino-induced muons is reversed relative 
to that for atmospheric muons. There are more 
negative than positive neutrino-induced muons 
{fj,~^ / fj,^ <1 / 2) , whereas for muons produced in 
the atmosphere, fi^ / fi'' ~ 1.25 This is a con- 
sequence of the ratio of cross sections, o^jay > 1, 



and the fact that — > /U while i^^ — > /Lt^ 
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Table 1 

Event rates for vertex contained and upward ex- 
ternal neutrino induced muons at the location of 
MINOS for the epoch of solar minimum in four 
energy intervals, with and without oscillations. 
The oscillation parameters used are sin^ 29 = 1 
and Am^ = 0.0025 eV^. The units for vertex con- 
tained events are lQ~^^g~^ .s~^ and for external 
upward going; muons arc W^^'^cm^'^ .s^^ 



E, GeV 


no oscillations 




0.0025 eV^ 




M 












contained vertex 


1 - 5 


6.88 


3.27 


5.14 


2.45 


5 - 10 


0.86 


0.41 


0.66 


0.31 


10- 20 


0.48 


0.24 


0.38 


0.18 


>20 


0.49 


0.23 


0.46 


0.21 






external 


upward 


going 


1 - 5 


3.57 


1.71 


2.01 


0.94 


5 - 10 


1.58 


0.75 


1.05 


0.48 


10 - 20 


1.55 


0.72 


1.21 


0.54 


> 20 


5.38 


2.19 


5.20 


2.10 



addition, the ratio y^/v^j, increases slowly above 
its low energy value of unity, first as muon-dccay 
becomes unimportant and at higher energy be- 
cause of the increasing importance of the channel 

p AK+ 

Secondly, the energy spectrum of the external 
muons is significantly harder than for the con- 
tained vertices. This is a consequence of the in- 
crease of the muon range with energy. A much 
higher range of neutrino energies contributes to 
the external events than to the events with in- 
teraction vertex inside the detector. For Ei, > 
100 GeV kaons become the dominant source of 
neutrinos, with a big contribution from the pro- 
cess p AK^ i^ij,- As a result, the charge 
asymmetry is somewhat larger for external events 
than for contained vertices. 

It is also interesting to look at the angular de- 
pendence of neutrino-induced muons with and 
without oscillations. Fig. 1 shows the zenith 
angle dependence of the neutrinos that produce 
muons with contained vertices and upward going 
neutrino-induced muons originating outside the 
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Figure 1 . a) Zenith angle dependence of the neu- 
trinos generating vertex contained muons of en- 
ergy above 1 GeV (upper line) and 10 GeV (lower 
line) and of upward going external muons (his- 
togram) in the absence of oscillations. The units 
for vertex contained events are {g.s.sr)~^ and for 
external upward going muons are (cm? .s.sr)~^ . 
The flux of upward going muons is divided by 10**. 
Note that the upward going muons are plotted up 
to cos 6 = 0.1; b) Same in the presence of oscil- 
lations with sin^ 26^ = 1 and Am^ = 0.0025. The 
dashed lines show the zenith angle dependence of 
the muons. 



detector. The strongest effect of oscillations is on 
the contained vertex events above 1 GeV. 

The symmetric shape of the no oscillation case 
is replaced by an asymmetric distribution with a 
deficit due to oscillations for upward muons from 
neutrinos with pathlengths comparable to the ra- 
dius of the Earth. This effect is in reality ob- 
scured by the large angle between the neutrino 
and muon direction at low energy as shown by the 
dashed lines in Fig. lb. The effect of the oscilla- 
tions is also obvious in the higher energy (E > 10 
GeV) sample. The rate of vertically upgoing con- 
tained vertex muons is decreased by a factor of 
^ 2. The muon-neutrino angle in this case is sig- 
nificantly smaller and does not affect much the 
oscillation features. The smallest effect is in the 
case of external events. On one hand, the contri- 
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Figure 2. Charge ratio of neutrino-induced muons 
a) for events with contained vertices; b) for ex- 
ternal muons. Solid lines are with no oscillations 
and dotted lines with oscillations for three bins of 
zenith angle centered on cos — —0.95, —0.55 and 
—0.15. Near vertical is highest solid line (right- 
most dotted peak). 



bution of high energy, non-oscillating neutrinos is 
the largest and on the other the downgoing ex- 
ternal muon rate is dominated by atmospheric 
muons, so most of the downward hemisphere is 
not available as a probe of neutrino oscillations 
with external events. The muon-neutrino angle 
for this sample is negligible. The background of 
atmospheric muons is low enough at the depth of 
Soudan so that neutrino-induced, external muons 
can be measured up to cos w 0.1, about 6° above 
the horizon. Integrated over all upward-going ex- 
ternal events the, oscillations decrease the rate 
by 22%, which is comparable to the uncertainty 
in the expected flux. [|l3| Instead, the effect of 
oscillations is visible as a distortion of the angu- 
lar distribution, as observed at MACRO |l^ and 
Super-K. || 

We assume here that neutrinos and anti- 



neutrinos have identical oscillations properties. 
Because of the properties of the differential 
charged current cross sections of neutrinos, how- 
ever, the relation between E^^ and is differ- 
ent for and u^. As a consequence, when 
the muon charge ratio is plotted as a function 
of muon energy for a given angular bin (corre- 
sponding to a given pathlength) oscillation fea- 
tures appear in the ratio because the oscillation 
minima show up at different energies for /i"*" and 
/i^. Figure 2 illustrates this for contained ver- 
tices (a) and for external muons (b). The ratios 
are shown for three bins of cos 9 going from nearly 
vertical (—1 < cos6 < —0.9) to nearly horizontal 
(-0.2 < cos6' < -0.1). The solid lines show the 
charge ratios in the absence of oscillations. The 
baseline charge ratio increases slightly at high en- 
ergy and is larger near the vertical because of 
the increased importance of kaon decay relative 
to pion decay. 

Calculation of external events involves an inte- 
gral over all contributing neutrino energies taking 
account of the increase of muon range with en- 
ergy. The growth of the charge ratio with muon 
energy is stronger than for contained vertices be- 
cause of the contribution of higher energy neu- 
trinos. In the oscillation case the variations are 
smaller than for contained vertex events and do 
not show the secondary sinusoidal features. For 
the external events, the peaks are wider and shal- 
lower than for contained vertices because of the 
contribution of a broader range of neutrino ener- 
gies. 

For more insight into the origin of this be- 
havior, we show in Fig. 3 the charge ratio with 
and without oscillations for the vertically upward 
bin of zenith angle, for which the pathlength 
L ^ 10"* km. The lower panel shows the corre- 
sponding survival probability, which is assumed 
identical for neutrinos and antineutrinos. An ap- 
proximate, quantitative understanding of the fea- 
tures can be obtained by considering the leading 
term in the differential neutrino charged current 
cross sections. For neutrinos 
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Figure 3. Upper part shows the variation of the 
fi~ / jj,'^ charge ratio of almost vertical {cos 9 ~ 
-0.95), vertex contained upward-going neutrino- 
induced muons as a function of muon energy 
for no oscillations (solid) and oscillations (dotted 
line). The lower part shows the neutrino survival 
probability for the same angular bin. 
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for antineutrinos. Here < z — Ef^/Ei, < 1 and 
CTo ~ 0.8 X lO^^^cm^. Using the approximation |l| 
we can estimate the relation for neutrinos by eval- 
uating the average E^ for the distribution 
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is a power-law approximation for the differential 
spectrum of neutrinos. Then 
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The corresponding result for antineutrinos is 

7 + 2 



7+1 



E^+. 



(6) 



With 7 « 2 in the relevant energy range Q , we 
have Ep/E^+ « 4/3 and E^/E^- « 2. Thus the 
first oscillation minimum at « 24 GeV should be 
reflected at w 18 GeV in and /z"*" as an increase 
in the ratio and at w 12 GeV in the op- 

posite channel, reflected as the minimum in the 
/Lt~//U+ ratio below 12 GeV. The slight deviations 
from these numbers in Fig. 3 reflect the breadth 
and asymmetry of the oscillation features folded 
with the neutrino spectrum. 

Unfortunately, the MINOS detector is not large 
enough to see the oscillation effect in the charge 
ratio well, if at all. For example, we estimate that 
there will be only some 12 /i~ and 5 /i"*" events 
with energy between 10 and 30 GeV in 25 kt-yr in 
the angular region — 1 < cos 9 < —0.5, where the 
charge ratio is high (see Fig. 2). Our estimates 
are based on the angular distribution of the in- 
teracting neutrinos and do not account for the 
experimental efficiency. 
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Figure 4. Charge ratio of ii^ / fi'^ for upward 
going contained vertex neutrino induced muons 
as a function of L/E for different zenith angles. 
Different symbols show events in cos 9 bins of - 
0.15, -0.25, -0.35, ... , -0.95. The horizontal case 
is not plotted because of the very strong change 
of the oscillation length which requires a much 
more exact calculation. 
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The best way to use all the data is to bin ev- 
ery event by L/E. Fig. 4 shows the charge ratio 
in upward going contained vertex events for the 
standard oscillation parameters as a function of 
L/E for nine different cosO bins. The horizontal 
bin is not included because the fast change of the 
oscillation length there requires a much more ex- 
act calculation. One can easily identify the peaks 
and valleys in the charge ratio that follow the 
oscillation structure of the atmospheric neutrino 
flux. The main feature is at L/i? of ~ 600, where 
the charge ratio can grow to almost twice its base- 
line value of yiT I [if^ « 2. The valley is at LjE 
of ~ 1200. One can subdivide the LjE range in 
two parts with different charge ratios, 2.36 in the 
LjE from 300 to 900 and 1.84 from 900 to 1500. 
The expected number of muons for 25 kT.yr are 
26.7 [T & 11.3 /Lt+ and 15.5 & 8.4 respectively, 
still very small statistics. 

One should also note that Fig. 4 shows the neu- 
trino pathlength L which is not experimentally 
measured. For some of these events, where the 
sinusoidal variation happens at energies below 10 
GeV, the angular smearing of the neutrino-muon 
angle will seriously obscure the effect when plot- 
ted as L^/i?^. 

It is interesting to ask if one could distinguish 
oscillations from neutrino decay with an obser- 
vation of the muon charge ratios. The LjE < 
1500 part of Fig. 4 would not change for a decay 
scenario that fits the muon disappearance obser- 
vations. At higher L/E values the charge ratio 
behavior would be much smoother without the 
secondary and tertiary peaks at L/E of about 
2000 and 3000. The increase of the charge ra- 
tio at these positions are however small and even 
more difficult to detect. 

The reliable detection of the variation of the 
muon charge ratio would require a much bigger 
detector. The MONOLITH (l6) detector with 
mass of 30 kT was proposed for the Gran Sasso 
Laboratory. If this detector had been built and 
were operated for five years it would be able to 
measure the effect at the 3cr level. 

The MINOS detector will be able to measure 
the total neutrino energy but assume that the 
muon and neutrino directions coincide. The dif- 
ference between the neutrino and muon energy 



spectra as a function of L, as well as the charge ra- 
tios, could provide additional observable param- 
eters for the more exact derivation of the oscilla- 
tion parameters. 

The exact manifestation of the effect however 
has systematic uncertainties related to the uncer- 
tainties of the atmospheric flux. These arise from 
the limited knowledge of the primary cosmic ray 
spectrum and the production of pions and kaons. 
The overall uncertainty for Ei, ^ 10 GeV is esti- 
mated in Ref. as ±25%. 

The uncertainty in the ratio i^p./v^^ should be 
significantly less than the uncertainty in the mag- 
nitude of the flux because the normalization of 
the primary cosmic ray spectrum cancels. The 
largest remaining source of uncertainty in the ra- 
tio is from production of charged pions and kaons, 
especially in the forward fragmentation region, 
as reflected by the spectrum weighted moments, 
Zp^+ , Zp^- , Zpx+ and Zp^^- . For pions and espe- 
cially for kaons, the positive channel dominates, 
and kaons become relatively more important at 
high energy. 

We can estimate the uncertainty in the i^p,/i^p, 
ratio by combining in quadrature the estimates in 
Ref. of uncertainty in the flux of + Vp from 
each of the four channels independently. This 
leads to an estimate of ±8% below 10 GeV and 
±9% for 10 < < 100 GeV. Analysis of data 
from forthcoming hadron production experiments 
with CERN |l|,|l| and Fermilab [|o) should lead 
to a reduction of this estimate. 

The effect we described above should also be 
detectable in a slightly different form in the 
long baseline oscillation experiment with the 
NuMI neutrino beam. The experimental statis- 
tics would be fully sufficient and the exact mani- 
festation would strongly depend on the beam en- 
ergy spectrum, purity, and composition. All of 
these three parameters are currently still at the 
design stage. 
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